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Abstract 

We derive extremal black hole solutions for a variety of four dimensional models 
which, after Kaluza-Klein reduction, admit a description in terms of 3D gravity coupled 
to a sigma model with symmetric target space. The solutions are in correspondence 
with certain nilpotent generators of the isometry group. In particular, we provide the 
exact solution for a non-BPS black hole with generic charges and asymptotic moduli 
in M = 2 supergravity coupled to one vector multiplet. Multi-centered solutions can 
also be generated with this technique. It is shown that the non-supersymmetric solu- 
tions lack the intricate moduli space of bound configurations that are typical of the 
supersymmetric case. 
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1 Introduction 



Soon after the attractor mechanism was first discovered in supersymmetric (BPS) black holes 
PQ, it was reformulated in terms of motion on an effective potential for the moduli [2]. Fer- 
rara et al demonstrated that the critical points of this potential correspond to the attractor 
values of the moduli. More recently, several groups used the effective potential to show that 
non-supersymmetric (non-BPS) extremal black holes can also exhibit the attractor mecha- 
nism, thereby creating a new and exciting field of research [31 H] . Many connections between 
non-BPS attractors and other active areas of string theory soon revealed themselves. An- 
drianopoli et al found that both BPS and non-BPS black holes embedded in a supergravity 
with a symmetric moduli space can be studied using the same formalism, and they uncov- 
ered many intricate relations between the two (5j[6]. Dabholkar, Sen and Trivedi proposed 
a microstate counting for non-BPS black holes (albeit subject to certain constraints [7j). 
Saraikin and Vafa suggested that a new extension of topological string theory generalizes 
the Ooguri-Strominger-Vafa (OSV) formula such that it is also valid for non-supersymmetric 
black holes [8]. Studying non-BPS attractors could also give insight into non-supersymmetric 
flux vacua. Given all these possible applications, it is important to characterize non-BPS 
black holes as fully as possible. 

There has been a great deal of progress in understanding the near-horizon region of these 
non-BPS attractors. The second derivative of the effective potential at the critical point 
determines whether the black hole is an attractor, and the location of the critical point 
yields the values of the moduli at the horizon; in this way, one can compute the stability 
and attractor moduli for all models with cubic prepotential [HI [TO]. However, the effective 
potential has only been formulated for the leading-order terms in the supergravity lagrangian. 
If one wants to include higher- derivative corrections, one can instead use Sen's entropy 
formalism, which incorporates Wald's formula, to characterize the near-horizon geometry in 
greater generality [TT] . Sen's method has led to many new results [T21 E31 Ell E51 ESI EZ|- 
The tradeoff is that this method cannot be used to determine any properties of the solution 
away from the horizon. 

The BPS attractor flow is constructed from the attractor value z* BPS = z^p^p 1 ,qi) by 
simply replacing the D-brane charges with the corresponding harmonic functions: 

zbps(x) = zgp^p 1 -> H : {x), q 1 -> #/(£)) (1.1) 

where the harmonic functions are 

/ H\x) \ = (h'\ ^fp 1 } 
\ Hjix) J \hj ) + \x\\ qi ) 

Moreover, this procedure can be applied to construct the multi-centered BPS attractor flow 
that describes the supersymmetric black hole bound state [18J, where the harmonic functions 
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are generalized to have multiple centers: 

/ H\x) \_fh I \ v 1 / {p J ) t \ 

It is conjectured in [19] that the non-BPS flow can be generated in the same fashion, namely, 
by replacing the charges in the attractor value with the corresponding harmonic functions. 
However, as will be proven in this paper, this procedure does not work for systems with 
generic charge and asymptotic modulijj 

In principle, one could construct the full non-BPS flow (the black hole metric, together 
with the attractor flow of the moduli) by solving the equation of motion derived from the 
lagrangian. However, this is a second-order differential equation and only reduces to a first- 
order equation upon demanding the preservation of supersymmetry. Ceresole et al have 
written down an equivalent first-order equation in terms of a "fake superpotential" , but 
so far, the fake superpotential can only be explicitly constructed for special charges and 
asymptotic moduli [211 [20]. The most generic non-BPS equation of motion is complicated 
enough that it has not yet been solved. Similarly, multi-centered non-BPS black holes have 
not been studied. 

Our goal is to construct the full flow for non-BPS stationary black holes in four di- 
mensions. Instead of directly solving the equation of motion, we reduce the action on the 
timelike isometry and dualize all 4D vectors to scalars. The new moduli space M.zd contains 
isometries corresponding to all the charges of the black hole, and the black hole solutions 
are simply geodesies on M.^d- This method was introduced in [22] and has been used to 
construct static and rotating black holes in heterotic string theory [231 121] and to study the 
classical BPS single-centered flow and its radial quantization [25], [26] . 

In this paper, we work in two specific theories of gravity, but we expect that this method 
can be used for any model whose M.zd is symmetric. The basic technique is reviewed in 
more detail in Section 2. In Section 3, we show how this method works in a simple case: the 
toroidal compactification of D-dimensional pure gravity. Section 4 serves as an introduction 
to single-centered attractor flow in Af = 2 supergravity coupled to one vector multiplet, and 
Sections 5 and 6 are dedicated to constructing the full flows for both BPS and non-BPS 
single-centered black holes with generic charges. We find that they are generated by the 
action of different classes of nilpotent elements in the coset algebra. Both types of flows are 
shown to reach the correct attractor values at the horizon. In Section 7, the procedure is 
generalized straightforwardly to construct both BPS and non-BPS multi-centered solutions. 
We use a metric ansatz with a flat spatial slice and we are able to recover the BPS bound 
states described by Bates and Denef. Using the same ansatz, we are able to build non-BPS 
multi-centered solutions. Unfortunately, solutions generated this way turn out to always have 

1 This has also been shown in [2U] . 
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charges at each center which are mutually local. The last section reviews our conclusions 
and suggests possibilities for future work. 



2 Framework 

Here we outline the method we will use to construct black hole solutions. This method was 
first described in [22] . We first reduce a general gravity action from four dimensions down to 
three, and derive the equation of motion. We then specialize to certain theories which have 
a 3d description in terms of a symmetric coset space. In such situations, we can easily find 
solutions to the equation of motion. The solutions are geodesies (or generalizations thereof) 
on the 3d moduli space and they are generated by elements of the coset algebra. 



2.1 3D Moduli Space 

We will study stationary solutions in a theory with gravity coupled to scalar and vector 
matter. Let the scalars be z % and the vectors be A 1 . Then the most general ansatz for a 
stationary solution in four dimensions is: 

ds 2 = -e 2U {dt + u) 2 + e- 2U g ab dx a dx b (2.1) 
F 1 = dA 1 = d(A I {dt + u) + A 1 ) (2.2) 

where a,b = 1, 2, 3 label the spatial directions and bold font denotes three-dimensional fields 
and operators. Since none of the fields are time-dependent, we can compactify on the time 
isometry and reduce to three-dimensional space M.$d- This procedure is called the c*-map. 
In three dimensions, a vector is Hodge-dual to a scalar. The equations of motion for uj and 
the gauge fields allow us to define the dual scalars W and 0a 1 - 

We then obtain the 3d lagrangian in terms of only scalars 

£ = ^v/g(-- R + d a <p m d a rg mn ) (2.3) 
where <p n are the moduli fields 

0" = {^#,<^,^,0 A i} (2.4) 

and g ab is the space time metric and g mn is the metric of a manifold M.%d- The system is 
3D gravity minimally coupled to a nonlinear sigma model with moduli space M.%d- Next, 
we will find the equation of motion in this theory. 
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2.2 Attractor Flow Equation 



The equation of motion of 3D gravity is Einstein's equation: 

R ab - ig afc R = KT ab = K(d a( p m d b rgmn ~ ^aA^^Qn^ (2.5) 

and the equation of motion of the moduli is: 

V a V>™ + T n mp d a <p m d a <f = (2.6) 

For simplicity, we consider only the case where the 3D spatial slice is flat (it is guaranteed 
to be flat only for extremal single-centered black holes). Then the dynamics of the moduli 
are decoupled from that of the 3D gravity: 

R a6 = d a <p m d b r9mn = (2.7) 

In the multi-centered case, we need to solve the full equations for the moduli as functions of 
the 3d coordinates x. For single-centered solutions, the moduli only depend on r; to satisfy 
the above equations, the motion of the moduli must follow null geodesies inside M. 3 d- 

A generic null geodesic flows to the boundary of the moduli space M.^d- A single- 
centered attractor flow is defined as a null geodesic that terminates at a point on the U — > 
— oo boundary and in the interior region with respect to all other coordinates. This is 
guaranteed for the BPS attractor by the constraints imposed by the supersymmetry. To find 
the single-centered non-supersymmetric attractor flow, one needs to find a way to construct 
null geodesies and a constraint that pick out the ones that stop at this specific component of 
the boundary. In the next section, we will show that we can do this for models with special 
properties, and the method can be easily generalized to find the multi-centered attractor 
solution. 



2.3 Models with Symmetric Moduli Space 



The problem of finding such a constraint in a generic model is not easy. To simplify, we study 
any model whose moduli space is a symmetric homogeneous space: M. 3 d = G/H. When 
A4 3 d is a homogeneous space, the isometry group G acts transitively on Ai 3D . H denotes 
the isotropy group, which is the maximal compact subgroup of G when one compactifies 
on a spatial isometry down to (1,2) space, or the analytical continuation of the maximal 
compact subgroup of G when one compactifies on the time isometry down to (0, 3) space. 
The Lie algebra g has the Cartan decomposition: g = h © k where 

[h, h] = h [h, k] = k (2.8) 
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When G is semi-simple, the homogeneous space is symmetric, and 

[k, k] = h 



(2.9) 



The models with symmetric moduli space includes: /J-dimensional gravity toroidally 
compactified to four dimensions, certain models of 4D Af = 2 supergravity coupled to vector- 
multiplet, and all 4D Af > 2 extended supergravity. The entropy of the last two classes is 
U-duality invariant. In the present paper, we will only consider the first two classes, namely, 
the Z)-dimensional gravity toroidally compactified to four dimensions, and the 4D Af = 2 
supergravity coupled to riy vector- multiplet. 

The left-invariant current is 

J = M~ x dM = J k + J h (2.10) 

where M is the coset representative, and Jk is the projection of J onto the coset algebra k. 
The lagrangian density of the sigma- model with target space G/H is given by Jk as: 

L = Tr( J k A * 3 J k ) (2.11) 



The geodesic of the homogeneous space written in terms of the coset representative is 
simply 

M = M e kT/2 with kek (2.12) 

where M parameterizes the initial point, and the \ is for later convenience. A null geodesic 
has zero length: 

|A;| 2 = (2.13) 

Therefore, in a homogeneous space, we can find the null geodesies that end at an attractor 
point by imposing the appropriate constraint on the null elements of the coset algebra. 

Since M is defined up to the action of the isotropy group H, in order to read off the 
moduli fields from M in an H-independent way, we construct the symmetric matrix using 
the metric signature matrix So: 

S = MS M T (2.14) 

In all systems considered in the present paper, H is the maximal orthogonal subgroup of G 
with the correct signature: 

HS H T = S for VU G H (2.15) 

That is, the isotropy group H preserves the symmetric metric matrix So. Therefore, S is 
invariant under M — >■ MH with H e H. Moreover, as the isotropy group H acts transitively 
on the space of of matrices with a given signature, the space of possible S is the same as the 
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symmetric space G/H. That is, the moduli of G/H can be combined into the symmetric 
matrix S. And the current of S is 

J s = S^dS (2.16) 



It is easy to perform the projection onto the coset algebra k. The (generalized) orthogo- 
nality condition of the isotropy group H can be expressed in terms of the subalgebra element 
h which is in H = e h as 

hS + S h T = V/ieh (2.17) 

In other words, (hSo) is anti-symmetric: (hSo) T = — (h5o). Thus the coset algebra, being 
the compliment of h, can be defined as the k with (kSo) being symmetric: (kSo) T = (kSo), 
i.e. 

^^S^kSo \/kek (2.18) 

Therefore, the projection of an element g in g onto the coset algebra k is: 

gk= 9+ w (219) 

For the left-invariant current J = M _1 dM, the projection onto k is: 

J k = J + S ^ jTS o 1 (2 . 20) 



It is straightforward to show that the current constructed from S is related to the pro- 
jected left-invariant current Jk by: 

J s = S^dS = 2(S M T y 1 J k {S M T ) (2.21) 

The lagrangian in terms of S is thus L = ^Tr(Jg A * 3 J 5 ). That is, the lagrangian density is 

£ = Itt^-VS • S^VS) (2.22) 

which is invariant under the action of the isometry group G: 

S -> G^SG where GeG (2.23) 
and whose conserved current is: 

J = S-WS (2.24) 

where we have dropped the subscript S in J s , since we will only be dealing with this current 
from now on. The equation of motion is the conservation of the current: 

V • J = V • (S^VS) = (2.25) 
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We now specialize to the single-centered solutions: they correspond to geodesies in the 
coset manifold. The spherical symmetry allows the 3d metric to be parameterized as 

ds 2 3 = C(r) 2 dx 2 (2.26) 

Then the equations of motion involve the operator d r r 2 C(r)d r , and reduce to geodesic equa- 
tions in terms of a parameter r such that 

| = *7<r). (2.27) 

The function C(r) is then determined from the equations of motion of 3d gravity. The 
equations of motion can be written as 

=0 (2.28) 
In the extremal limit the geodesies become null, the 3d metric is flat and 

t = — (2.29) 

r 



In the search for multi-centered extremal solutions, where the spherical symmetry is 
absent, it is very convenient to restrict to solutions with a flat 3d metric. This is consistent 
with the equations of motion as long as the 3d energy momentum tensor is zero everywhere: 

T ab = Tr(J a J h ) = (2.30) 

The coupled problem with generic non-flat 3d metric is much harder, and exact solutions 
are hard to find unless a second Killing vector is present. 



Since different values of the scalars at infinity are easily obtained by a G transformation, 
to start with, we will consider the flow starting from Mq = 1, and generalize to generic 
asymptotic moduli later. For a single-centered solution, the flow of M is M = M e kr ^ 2 . 
Since all the coset representatives can be brought into the form e 9 with some g G k by an 
H-action, we can write M = e s/2 , so M = e 9/2 e kr/2 . And the flow of S is 

S(t) = e 9 / 2 e kT e 9 ' 2 S Q (2.31) 

The charges of the solution are read from the conserved currents 

;(r) = f7^^^ (2.32) 
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3 Torus Reduction of D- dimensional Pure Gravity 



Now we use the method introduced in the previous section to analyze pure gravity toroidally 
compactified down to four dimensions. We explain why the attractor flow generator, k, 
needs to be nilpotent, and we find the Jordan forms of k 2 and k. Using this information, we 
construct single-centered attractor flows. We then generalize to multicentered black holes 
in pure gravity and show that these solutions have mutually local charges and no intrinsic 
angular momentum. 



3.1 Kaluza-Klein Reduction 



The simplest example of a system that admits a 3d description in terms of a sigma model 
on a symmetric space is pure gravity in D dimensions, compactified on a D — 4 torus. The 
KK reduction to AD parameterizes the metric as 

ds 2 D = Ppq (dy p + A*daf)(dy q + A%dx v ) + ^=^4 1 < p, q < D - 4 (3.1) 

Here y p are the torus coordinates, x^ the coordinates on R 3,1 , and p pq the metric of the torus. 
A 4D metric with one timelike Killing spinor is then parameterized as 

ds\ = -u(dt + Uidx 1 ) 2 + ^ds 2 3 (3.2) 

where u = e 2U , to connect with the parametrization in the later part of the paper; and 
i = 1, 2, 3 denote the 3d space coordinates. 

The two expressions combine as 

ds 2 D = G ab (dy a + ufdx*) (dy b + tfdx^ + — 7—pzdsl 0<a,b<D-4 (3.3) 

J — det G 

Here y a are the torus coordinates plus time, x % coordinates on R 3 and 

r — ( Ppq PprA r \ , s 

and uo a = (Cu p , uj°) is: 

Co p = (A p - A p u,i)dx l Co° = cu (3.5) 

If the forms u a are dualized to scalars a a as 

da a = — det GG ah * 3 duo b (3.6) 



10 



the various scalars can be combined into a symmetric unimodular (D — 2) x (D — 2) matrix 



S = 

In terms of the 4D fields that is 



S 



Ppq - 



1 



V 



Uy/det p 

1 --a 



uy^detp 1 



G a b + 



1 



detG 



a a a b 



i 



detG 



«6 



dotG Qa 
1 

detG 



a a q A r p rs A s - 



i 



ity'detp 



a p a 



i 



v / dctp 



u^fdetp 



a a 



it-^dct p P 
U\/dct p^® 
u\/det p ) 



(3.7) 



(3. 



The equations of motion derive from the lagrangian density £ = TrV S ^ _1 V 5 S^ 1 , in- 
variant under S — > U T SU for any U in SL(D — 2). As this SL(D — 2) action is transitive 
on the space of matrices with a given signature, the space of possible S is the same as 
the symmetric space SL(D — 2)/SO(D — 4, 2). Notice that the signature of the stabilizer 
SO(D — 4, 2) is appropriate for the reduction from (D — 1,1) to (3, 0) signature. The usual 
reduction from (D — 1, 1) to (2, 1) would give a SL(D — 2)/SO(D — 2), while the Euclidean 
reduction from (D, 0) to (3, 0) gives SL(D - 2)/SO(D - 3, 1) [27]. 



The coset representative under the left SO(D — 4, 2) action can be described in terms of 
a set of vielbeins 

e A = E*{dy a + ufdx 1 ) e 1 = — ^ef^ (3.9) 



as 



M 



1 -,oc h 1 



det_E"° det_E 

Then the symmetric SO(D — 4, 2) invariant matrix 

S = MS M T 



(3.10) 



(3.11) 



can be used to read off the solution more easily. Without loss of generality, we can take So 
to be the signature matrix: 



S = Diag(r), -1) = Diag(l, ■■■ ,1,-1, -1) 



(3.12) 



The equations of motion are equivalent to the conservation of the SL(D — 2) currents 
J = S~ 1 dS. Some of those currents correspond to the usual gauge currents in 4D: the first 
D — 4 elements of the last column are the KK monopole currents, the first D — 4 

elements of the row before the last Jd-3,i are the KK momentum currents and the element 
Jd-3,d-4 is the current for the 3d gauge field u. Regular AD solutions must have zero sources 
for this current, otherwise uj will not be single valued. 
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3.2 Nilpotency 



We now show that all the attractor flows are generated by the nilpotent generators in the 
coset algebra. To get extremal black hole solutions with a near horizon AdS 2 x S* 2 , the 
function u must scale as r 2 as r goes to zero while the scalars go to a constant. This makes S 
diverge as ^. The most natural way for S to diverge as r 2 for large r is that k is nilpotent, 
with 

k 3 = . (3.13) 
This will be the crucial condition through the whole paper. 



3.3 A Toy Example: Hyperkahler Euclidean Metrics in Four Di- 
mensions 



Not every null geodesic corresponds to extremal black hole solutions. Let's consider a simple 
example: hyperkahler euclidean metrics in AD. 

Although this example is not about a black hole, it is still quite instructive. The 3d sigma 
model is SL(2)/SO(l, 1), i.e. AdS 2 - The coset representative is written as 

M = ( I 2 1) (3,4) 

V u 1 / 2 u l/2 J 



and the symmetric invariant 



u u 



S=[ U a « Y ) (3-15) 



A geodesic is the exponential of a Lie algebra element in the orthogonal to the stabilizer. 
The stabilizer 5*0(1, 1) is generated by a 1 . A null geodesic is hence the exponential of 
k = a 3 ± io 2 . This is a nilpotent matrix, k 2 = 0, hence M — 1 + rk/2. Take: 

k = a 3 + ia 2 = ^ _^ _^ ^ (3.16) 

then 

and we can read off the geodesic solution from the invariant 
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Hence 

«=—!—, a = --^—. (3.19) 

1 — r 1 — r 

Dualizing a, we get 

da 1 , i _ 1 .„ 

* duo = = — dr u^ 1 = 1 + - (3.20) 

u r 2 r 

The 4D Euclidean metric is just the Taub-NUT metric. Notice that the multi-centered 
Taub-NUT generalization of the metric is obtained by replacing r above with some harmonic 
function £V r^i ■ The sigma model equations of motion are equivalent to the conservation 

of the current J = S' _1 VS', and if 5 is given as above with r = t(x) then the equation of 
motion are 

V 2 r(f) = (3.21) 



3.4 Single-centered Black Holes in Pure Gravity 
3.4.1 Constructing the flow generator k 

Now we look in detail at the single-centered black holes in pure gravity. Notice that as u 
goes to zero S tends asymptotically to a rank one matrix 



1 



U\/det p 



ao<y q aoao ao ) (3.22) 
a q ao 1 



hence the matrix k 2 should also have rank 1. By inspection of S it is clear that a k 2 of rank 
higher than one gives a geodesic for which the matrix elements of p also diverge as r 2 so that 
the scalar fields do not converge to fixed attractor values. 

Notice that if k is nilpotent, then S is a polynomial in and the various scalars in the 
solution will all be simple functions of r for such extremal solutions! 

The explicit form for k in terms of the charges is then straightforward to write. Consider 
the Jordan form of k 2 : as it is nilpotent, the eigenvalues are all zero. As it is rank one, it 

1 



has one single indecomposable block of size two: ^ ^ ^ J . It is written as 

k 2 = -r]vv T f(P,Q) (3.23) 

with v null in the metric rj, and f(P, Q) any degree- two homogeneous function of the charge 
(P, Q). This form is chosen so that v does not scale with the charge (P, Q). 
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Then k must have a Jordan form with all eigenvalues zero, one block of size 3:1001 

\ ) 

and possibly some other extra blocks of size two. Alternatively, there is a subspace V anni- 
hilated by A;, a subspace V whose image under k sits in V and has the same dimension as 
V, and a single vector w such that kw G V and is non-zero. 

kw c V, W c V, kV = 0. (3.24) 

From the symmetry of r]K, it follows that the space kV is made of null vectors only, and 
that kw is orthogonal to it. Because r/k 2 = —vv T , (kw) T rjkw is negative. Because of the 
signature of rj it is straightforward to see that V can be of dimension at most one; hence 
there are no blocks of size 2 in the Jordan form of k. 



Taking this into account, the final form of k is simply 

k = r]wv T + 7]vw T (3.25) 

where v and w are two orthogonal (D — 2)-dimensional vectors with v being null and w 
having norm —f(P, Q)'- 

wr/v = 0, v rjv = 0, WTjw = —f(P,Q). (3.26) 

Using the fact that in k, the first D — 4 elements of the last column K it D_2 are the magnetic 
charges, and the first D — 4 elements of the row before the last K£,_ 3 i are the electric 
charges, and the element -Kd-3,d-4 is the Taub-NUT charge, which has to vanish, we have 
2(D — 4) + 1 = ID — 7 conditions. Together with the three constraints coming from the 
norms and orthogonality condition f)3.26p . they can be used to solve for the 2D — 4 degrees 
of freedom in (v, w). 



The full solution of (v, w) requires one to solve some degree-four equations, hence we'll 
leave it in a slightly implicit form. Let (p, q) be two (D — 4)-dimensional vectors proportional 
to the magnetic and electric charges, so that the magnetic charge and the electric charge 
(P, Q) of the 4D gauge fields are 

P = \Jp 2 + p ■ qp Q = \/q 2 + p ■ qq (3.27) 

And we choose f(P, Q) to be 

f(P,Q)=p.q (3.28) 
The solution of v and w written in terms of (p, q) is 

/ q + p \ / —— I q-p 

f -^qi+p.q , w= Vl-A\ -vV+p-g I. (3.29) 



^ P ' \ -\/p 2 + p- q ) \ \/p 2 + p ■ q 
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and k can be written as 



T T 

qq • — pp 
y/ q 2 +p-~qq T 
—\/p 2 +p- qp T 



-yV +p - qq \Jp 2 + p- qp 
-(q 2 +p-q) 

p 2 + p ■ q 



(3.30) 



3.4.2 Full flow 



First, for the full flow starting from M = 1, the scalars for the attractor solution generated 
by this k can be read off from S(t) = e kT , by comparing with the form of S in terms of the 
4D fields: 

u- 2 = [l + (p 2 + p.q)( T + ^r 2 )][l + (q 2 + p • q)(r + ^-r 2 )} (3.31) 

and 

(gg T - pp T )r + {(p 2 + p ■ q)qg T - g(p + q)(p + g) T ]r 2 
P ~ + l + (p 2 + p-q)(r + ^) {6 - 6Z) 

Notice that as r — > oo, T~ 2 e~ 2U has the correct limit which is the entropy where P and 
Q are the physical electric and magnetic charges. 

To generalize to arbitrary asymptotic moduli, M(r) = e 9 ^ 2 e kr ^ 2 , and the flow of S is 
fl2.3ip . which can be written as S(r) = 6 K ^Sq, where K(r) is a matrix function. From now 
on, we use lower case k to denote the coset algebra that generates the attractor flow, and 
capital K to denote the function which we exponentiate directly to produce the solution. 

We will choose K(r) to have the same properties as the generator k: 

K 3 (t) = and K 2 (t) rank one (3.33) 

The equations of motion V • (S^VS*) = then simplify considerably with this ansatz. If 
one further requires that the subspace image of K 2 (t) remains constant everywhere, such 
that 

K 2 (t)VK(t) = VK(t)K 2 (t) = (3.34) 

then the current reduces to 

J = S^VS = S ^VK(r) + ^[VK(t), K(t)}^J So (3.35) 

and the equations of motion are 

V 2 K{t) + ^[V 2 K(t), K{t)} = (3.36) 
which is solved by a harmonic K{r). 
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It might appear hard to find a K(t) that is harmonic and satisfies all the required 
constraints. However, by remembering that the constraints dictate K{r) to have the form: 

K(t) = r]VW T + 7]WV T (3.37) 

with V being null and doesn't scale with the charge (P,Q), and W orthogonal to V ev- 
erywhere, one can simply pick a constant null vector V — v' and a harmonic vector W(r) 
everywhere orthogonal to v'\ 

W(t) =w'r + m with v ■ W(r) = (3.38) 

Here m is a (D — 2)-vector and contains the information of asymptotic moduli. Thus an 
appropriate K{r) is built: 

K(t) = k'r + g (3.39) 

where 

k! = r]v'w' T + 7]w'v' T g = rjv'm T + r\mv T (3.40) 

Now we need to solve for (v',w') for the same charge (P,Q) but in the presence of m. 
The form of g guaranteed that 

[k',g] = (3.41) 

where we have used the fact that v' is null and w' is orthogonal to v'. Therefore, shifting the 
starting point of moduli does not change the current as a function of (v,w): 

jw, = s„ g) s„ = so { n^r^r ^ So (3 . 42) 

Thus, the solution of (v 1 , w') in terms of charges solved from the current does not change as 
we vary the starting point of the flow, i.e. they do not depend on the asymptotic moduli: 

v'(Q) = v(Q) w'(Q) = w(Q) (3.43) 

In summary, the flow with arbitrary starting point is simply generated by 

K(r) = 7]vW{t) t + t]W{t)v t with W = wr + m (3.44) 
where (v, w) only depend on the charges (P, Q) and m gives the asymptotic moduli. 



3.4.3 Example: 5D pure gravity compactified on a circle. 



Consider for example the case of extremal black holes in D = 5 pure gravity compactified 
on a circle. The 3d sigma model is SL(3)/SO(l, 2). The symmetric invariant is 



\ 



-^"i^o p(A ) 2 



u 

_± 



F ai \ 

1 



a 



(3.45) 
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Then we can calculate S = S e kT using (13.301) and compare the result to S gr above to solve 
for all the scalars. We find that 



e~ 2U = A /[l + ( g 2 + M)(r + ^ r2)][1 + (p2+M)(r + ^ r2)] (3 . 46) 



P = 1 I I I Sfe I p\J (3-47) 



l + (q 2 +pq)(T + fr^ 
1 + (p 2 + pq)(r + fr 2 ) 

when starting from the identity. If we allow arbitrary g the flow is too complicated to write 
explicitly here, but the attractor value of p is the same: q/p. 



3.5 Multi-centered Solutions in Pure Gravity 



In the context of pure gravity compactified on a torus, we can also give some examples of 
multi-centered solutions in the same spirit as the ones for BPS solutions in N = 2 super- 
gravity, though some important features of the latter are not present here. 



We are interested in solutions given in terms of harmonic functions which can generalize 
the single-centered extremal solutions presented above. Similar to the single-centered case, 
we exponentiate a matrix function K(x): 

S(x) = e K{S) S (3.48) 

We will choose K(x) to have the same properties of the generator k: 

K 3 (x) = and K 2 (x) rank one (3.49) 



Using a similar argument to the single-centered flow, we require that the subspace image 
of K 2 (x) remains constant everywhere, such that K 2 (x) VK(x) = VK(x)K 2 (x) = 0, then 
the equations of motion are 

V 2 K{x) + ^[V 2 K(x), K(x)} = (3.50) 
which is solved by a harmonic K(x). 

A multi-centered K{x) that is harmonic and satisfies all the required constraints can then 
be built in the same way as the single-centered one: 

K(x) = r]vW(x) T + rjW(x)v T (3.51) 

where v is the same constant null vector as in k, and W(x) is the multi-centered harmonic 
function: 

W(x) = V — ^— + m (3.52) 

i ' 1 ' 
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where x\ is the position of the zth center, and is determined by the charges at the ith 
center, and m is related to the moduli at infinity. Requiring W(x) to be orthogonal to v 
everywhere gives the following constraints on the {wi, m}: First, taking x to infinity, it gives 

vm = (3.53) 

Second, Wi are orthogonal to v: 

v-Wi = (3.54) 

In addition to the constraint from the zero Taub-NUT charge condition — bz — cy = 0, this 
makes the space of each possible Wj only (D — 4)-dimensional. That is, though W is a 
(D — 2)-vector, one has only (D — 4) independent harmonic functions to work with, because 
of the orthogonality to v and the requirement of no timelike NUT charges. This makes the 
solution relatively boring. 

The multi-centered solution in pure gravity does not have the characteristic features of 
the typical BPS multi-centered solution in M = 2 supergravity, where many centers with 
relatively non-local charges form bound states which carry intrinsic angular momentum. 

The basic reason is that when the ansatz K(x) = r]vW T (x) + r]W(x)v T is used, the 
second term of the conserved currents J = So(VK + ^[VK, K])Sq drops out. The first 
result is that the charges of the various centers in the solution can be read off directly from 
(v,Wi), and they do not depend on the positions, charges of the other centers. Thus, there 
is no constraint on the position of each center as in the M = 2 BPS multi-centered solution; 
centers can be moved around freely. 

Moreover, the condition of no timelike Taub-NUT charge is a linear constraint on the 
charges at each center which results in a static 4D solution, as *du = leads to u = 0. 
Therefore, no angular momentum is present. 



4 Attractor Flows in G 2 (2)/(SL(2,R) x 5L(2,R)) 

We now tackle a more complicated subject: Af = 2 supergravity coupled to one vector 
multiplet. First, we reduce the theory down to three dimensions and derive the metric 
for the resulting moduli space, which is the coset G 2{2) /(SL(2, R) x SX(2,R))0 We then 
discuss the Cartan and Iwasawa decompositions of the group G 2 (2), which we use to construct 
the coset algebra and translate the flow of coset representative into the flow of the moduli 
fields, respectively. We then specify the representation of G 2 (2) we will be working with, and 
describe the form of attractor flow generators in this representation. 

2 Other work on this coset space has appeared recently, including p51 f^l 150] , 
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4.1 The moduli space A4%d 



The 3d moduli space for J\f = 2, d = 4 supergravity coupled to n v vector multiplets is well- 
studied, for example in [2H E21 E3] • Some of the main results are compiled in the Appendix. 
Here we briefly review the essential points. 

The bosonic part of the action is: 

S = --L J d 4 xy^ R - 2g i] dz i A * 4 dz~ j - F 1 A G, (4.1) 

where 1 = 0, l...ny, and Gj = (ReAf)ijF J + (ImAf)u * F J . For a model endowed with a 
prepotential F(X), 

(ImF ■ X) I (lmF ■ X)j , . 

^"- F " + 2i X-imF.X (4 ' 2) 
where Fjj = didjF(X). We reduce to three dimensions, dualizing the vectors (a;, A 7 ) to the 
scalars (a,Bj), and renaming Aq as A 1 . The metric of M.% D is then 

ds 2 = dU-dU + ^e~ 4U (da + A^Bj - B^A 1 ) ■ (da + A'dfi, - B^A 1 ) + g fj (z, z)dz l ■ dz j 

+^e- 2U [(ImAf- 1 ) IJ (dB I + N IK dA K ) ■ (dBj + Jf JL dA L )\ (4.3) 

It is a para-quaternionic-Kahler manifold. Since the holonomy is reduced from SO (Any + 
4)) to Sp(2, M) x Sp(2riy + 2,IR), the vielbein has two indices (a, A) transforming under 
Sp(2, M) and Sp(2riy + 2,R), respectively. The para-quaternionic vielbein is the analytical 
continuation of the quaternionic vielbein computed in [34]. The explicit form is given in the 
appendix. 

For ny = 1, X 1 = (X^X 1 ). For our purpose, we choose the prepotential 

f(x) = -^r ( 4 - 4 ) 

The metric of A4l D with one-modulus is (14.31) with g z2 = and M and (ImN)^ 1 being 
N = l ~(2x - iy)(x + iy) 2 3x(x + iy) \ ImA f-i _ _ 1 ( 1 x 



3x(x + iy) —3{2x + iy) J y 3 \ x 3x 2 + y 2 



The isometries of the Ad^ D descend from the symmetries of the 4D system. The gauge 
symmetries in 4D give shifting isometries of Ai^ D , whose associated conserved charges are: 

qi dr = J A i = P A i - BiP^ tfdr = J Bl = P Bl + A l P a , kdr = J a = P a 

(4.5) 
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where the momenta {P a , P^i, -Pb 7 } are 

P a = ^e" 4C/ (da + A'cLB/- B^A 1 ) (4.6) 

P A x = e- 2U [(ImM)udA J + (ReAf) IJ (ImAf- 1 ) JK (dB K + {ReM) KL dA L )\ - BjP^.7) 
P Bl = e- 2U [{ImU- 1 ) IJ {dB J +{ReM)j K dA K )]+A 1 P a (4.8) 

Here r is the affine parameter defined as dr = — *3 sm6d8d<j). (p°,p 1 ,qi,qo) are the D6- 
D4-D2-D0 charges, and k the Taub-NUT charge. A non-zero k gives rise to closed time-like 
curves, so we will set k = from now on. 

Note that the time translational invariance in 4D gives rise to the conserved current 

J U = P U + 2a J a + A*J A i + B^bi (4.9) 
where Pjj = 2dU. The corresponding conserved charge is the ADM mass: 2M.ADudT = Jjj. 



4.2 Extracting the Coordinates from the Coset Elements 



The metric (14.31) for the case ny = 1 describes an eight- dimensional manifold with coordi- 
nates (p n = {u,x,y,a, A , A 1 , Bi, B }. This manifold is the coset space G 2 (2)/ (SL(2, M) x 
SL(2, M)). The root diagram for the Cartan decomposition of £2(2) is shown in Figured! The 
six roots that lie on the horizontal and vertical axes {L^, L\, L^, L^} are the six non-compact 
generators of the subgroup H = SL(2,R) h x SL(2,R) V : 



(4.10) 



and the two vertical columns of eight roots {a a A} are the basis of the subspace K. {a^, a 2 A} 
for each A is a spin-1/2 doublet under the horizontal SL(2, M): 



[L 



CLia 



-\a>iA 

\a 2 A 



[L 



h ■ 



aiA 
a 2 A 







[L 



h ■ 



aiA 
a 2 A 



°2A / v 2 

And {a a i, a a2 , a a3 , a a {\ for each a span a spin-3/2 representation of the vertical SL(2, 



a 2 A 




[L 



( a a i \ 

\ a a4 / 



/ — §Oq:1 \ 
— 2 a a2 

2 G a4 / 



V 



( a a i \ 

\ 0>ai J 



( \ 

3(l a i 

2a a2 

\ a «3 / 



( a al \ 

\ a aA / 



/ -a a2 \ 
— 2a Q 3 





/ 



All the commutators can be easily read off from the Root diagram (pQ), we will only write 
down the following ones which will be useful later. 



an, au 



--[a 12, ai3j 



[a 2 i, a 2 4\ 



•-[a 2 2, ^23] 



-AT, 



(4.11) 
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3-22 • 



-•-hi. 



a 2 3 • 



a 12 • 



a i3 • 



a 2 4 • a 14 • 

Figure 1: Root Diagram of Cartan Decomposition of £2(2) 



Being semisimple, the algebra of £2(2) has the Iwasawa decomposition g = h©a©n, where 
a is the maximal abelian subspace of k, and n is the nilpotent subspace of the positive root 
space S + of a. In Figure [2J we show the Iwasawa decomposition of G 2 (2)- The two Cartan 
generators in a are {u, y}, and {x, cr, A , A 1 , B 1; B } span a nilpotent subspace n: n 7 = 
for n G n. a and n together generate the solvable subgroup Solv of G, which act transitively 
on M.3D = G2(2)/ SL(2, K) x SL(2,M). In particular, y generates the rescaling of y, and 
{u, x, cr, A , A 1 , Bi, B } generates the translation of {U, x, a, A , A 1 , B\, B }. The moduli 
space M.3D can be parameterized by the solvable elements: 



E(0) 



3 !7u+(ln3/)y xx+A 1 A 1 +B I B 1 +ao- 



The origin of the moduli space 

a = A = A 1 = B t = B = 
correspond to = 1. 



x = 



y = u 



(4.12) 



(4.13) 



In Fig [2j the isometries are plotted according to their eigenvalues under the two Cartan 
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A • 



A • 



-•-a- 



x O 



-a-»- 



A 1 • 



A • 



Figure 2: Root Diagram of Isometry of M 3D = G 2 p)/(SL(2,R) x SL(2, 
{u, y, x, cr, A , A 1 , Bi, B } generates the solvable subgroup. 



generators u and y [31j. {u, y} are related to a a A byE|: 



U = ~o + fl24 ) ~ ( ai 3 + fl 22)] 



y = ot 3 ( ai1 + fl24 ) + ( fll 3 + a 22)] ( 4 - 15 ) 



The three generators {cr, u, cr} on the horizontal axis and {x, y, x} on the vertical axis form 
the horizontal and vertical SL(2, M), respectively. The vertical SL(2, M.) generate the duality 
invariance. Denote the two vertical columns of eight isometries as 



/ 61 


€11 


\ 


/ 


-A 


B 


\ 










-3A 1 


3Bi 




€23 


lis 






B x 


-A 1 




V 


€14 


/ 


V 


-B 


A 


/ 



{£iA, $2a} for each A span a spin- 1/2 representation of the horizontal SL(2, 
for each a span a spin-3/2 representation of the vertical SL(2,M). 



3 The matrix representation of u and y are 



and{£ Q i,£ a2 ,£ 



u = Diag[0^, -i 0,-|, i 0] 



y = 2>ia ff [l,-i,-|,-l,|,| f 0] 



(4.14) 
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Parameterizing the coset representative M as the solvable elements, the symmetric matrix 
S can be written in terms of the eight coordinates 0™, from the solvable elements £ 

S(0) = S(0)S'oE(0) T (4.16) 

The coordinates are read off from the symmetric matrix S. 

4.3 Nilpotency of the Attractor Flow Generator k. 

The near-horizon geometry of the AD attractor is AdS 2 x S 2 , i.e. 



*r as t — > oo (4.17) 



-1/ 



In terms of the variable w = e 2i7 

u — > — 3— r~ 2 as r — > oo (4.18) 

The solvable element is 

M = e Uu+ ~~u* u (4.19) 
As the flow goes to the near-horizon, u — > 

M(r) ~ -u " £ / 2 ~ / (4.20) 

where — £ is the lowest eigenvalue of u. That is, M(r) is a polynomial function of r. 

On the other hand, since the geodesic flow is generated by k via 

M(r) = M(0)e kr/2 (4.21) 

i.e., M(t) is an exponential function of r. To reconcile the two statements, k must be 
nilpotent: 

k e+1 = (4.22) 

That is, the element in k that generates the attractor flow is nilpotent. Moreover, by looking 
at the weights of the fundamental representation of G 2 , we see that £ = 2 

kl = (4.23) 

Moreover, the nilpotency of the attractor flow generators guarantees that it is null: 

k 3 = =^ (k 2 ) 2 = Tr(k 2 ) = (4.24) 
which means that k is null. 
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4.4 Properties of Attractor Flow 



The scalar moduli space is parameterized by a symmetric 7x7 matrix S which sits in G 2 (2) , 
i.e. preserves a non-degenerate three form Wijk such that r] is = WijkW s t u w mno e^ ktumno is a 
metric with signature (4, 3) normalized so that if = 1. To facilitate the comparison with the 
pure 5.D gravity case we decompose 7 as 3 © 3 + 1 of SX(3) and pick as non-zero components 
of w the 3 A 3 A 3, 3 A 3 A 3 and 3 <g> 3 ® 1 as 

w — dxi A c&r 2 A dz 3 + cV A cfo/ 2 A dy 3 \=dx a A rfy a A dz (4.25) 

v2 

The resulting expression for r\ is 

7] = dx a dy a - dz 2 (4.26) 

We know that k must be an element of G 2 (2), hence also of 50(4,3). As in the pure 
gravity case, we choose the representation such that 

S kS = k T S k 2 S = (k 2 ) T (4.27) 

In this base a G 2 (2) Lie algebra element is given as 

/ All £nj2kV k V2wn \ 
k = e i2 ^ lk w k -Af 2 -V2v l2 (4.28) 
V -V2v^ V2w j2 / 

Here A is a traceless 3x3 matrix. S is a symmetric element in G 2 (2) with signature 
{1, -1, -1, 1, -1, -1, 1}, i.e. S = MS M T with 

77! 

50 = 1 77! | =^(1,-1, -1,1, -1,-1,1) (4.29) 
1 

where 771 is the one for pure gravity. 

If the gauge field is turned off, then S is block diagonal 



S\ F=0 = S- 1 (4.30) 




where S gr is the same as the one for pure 5D gravity. Turning on a non-zero 5D vector field 
corresponds to a more general S: 

S = e k T(S\ F=0 )e^ (4.31) 



24 



with k\ a Cj2(2) Lie algebra matrix with w\ equal to the fifth component of the gauge field, 
v 2 equal to the time component of the gauge field and w 3 equal to the scalar dual to the 
three-dimensional part of the gauge field. 

In this representation, (x, y) can be extracted from symmetric matrix S via: 

^33 IT) ^33 (^J 



And u via: 



u = , (4.33) 
V / ^33(r)5 55 (r) - S 35 (ry 



The 4D gauge currents sit in J = S^VS, where J^iJsi) is again the electric (magnetic) 
current for the KK photon, J 32 the timelike NUT current, and .772(^51) the electric (magnetic) 
current for the reduction of the 5D gauge field. 

J 32 = -2J CT Ji 2 = V2J A o J 72 = -J A i J 51 = -V2J Bl J 3 i = V2J Bo (4.34) 
Moreover, 

J22 - J33 = 2Ju (4.35) 

We use Q to denote the charge matrix, where it relates to the D-brane charge {p ,^ 1 , qi, qo} 
and the vanishing NUT charge k by 

(Q3i,Q 5 i,Q72,Qi2) = (v / 2p°,-v / 2p 1 ,^i,v / 2?o) Q 32 = -2A; = (4.36) 



Since k is nilpotent: k 3 = 0, 

S = e kT S =(l + kr + lk 2 r 2 )S (4.37) 

The AdS 2 x S 2 near-horizon geometry of the 4D attractor dictates u = v ^ r~ 2 as r — > 00. 
Therefore, the flow generator k can be obtained by 

k 2 = 2V BH UuS\ u ^ )S (4.38) 

Computing k 2 using S constructed from the solvable elements S(0) shows that k 2 is of rank 
two, its Jordan form has two blocks of size 3. It can be written as 

k 2 = v « v bC ab S (4.39) 

a,6=l,2 
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with v a null and orthogonal to each other: v a ■ v b = v^SqVi, = 0, and c a b depends on the 
particular choice of k. Thus k can be expressed as: 

k= J2( V a W a+ W « V a) S (4-40) 
a=l,2 

where each w a is orthogonal to both v a \ w a ■ v b = 0, and w a satisfy 

w a -w b = c ah (4.41) 

Parallel to the pure gravity case, the single-centered attractor flow is constructed as 
Sir) = e K ^S , where we choose K(t) to have the same properties as the generator k: 

K 3 (t) = and K 2 (t) rank two (4.42) 

This determines K(r) = kr + g where 

k = ^2i v a w I + w aV a T }S and g = ^ K m I + m a.vl\So (4.43) 

o=l,2 a=l,2 

where the two 7-vectors m a 's are orthogonal to f a and contain the information of asymptotic 
moduli. Using [[k, g],g] = 0, the current is reduced to 




(4.44) 



from which we obtain v a and w a in terms of the charges and the asymptotic moduli. 



5 Flow Generators in the G 2{2) /{SL{2 1 R)x5L(2, R)) Model 

We now explicitly construct the generators of single-centered attractor flows. We start with 
the BPS flow which is associated with a specific combination of the coset algebra generators 
a a A- It can be derived from the condition of preservation of supersymmetry. We then 
construct the non-BPS attractor flow generator in analogy with the BPS one. In Section 
5.2, we write ksps an d k nonBPS in terms of the v a and w a vectors. This form will be especially 
helpful in generalizing to the multi-centered case. 



5.1 Construction of Flow Generators 
5.1.1 Constructing k BPS using supersymmetry 

To describe BPS trajectories it is useful to remember that the stabilizer of S in G 2 (2) is 
SO(l,2) x 5*0(1,2), corresponding to the elements of G 2 (2) which are antisymmetric after 
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multiplication by Sq. Geodesies are exponentials of elements that are symmetric after mul- 
tiplication by Sq. Such elements sit in a (2, 4) representation of 50(1, 2) x S0(1, 2). A BPS 
trajectory is highest weight for the first S0(1, 2). Labelling the symmetric generators as a aA 
under the two S0(l,2) groups, a BPS trajectory is generated by 

k B ps = a aA C A z a . (5.1) 

The twistor z and the coefficients C A are fixed in terms of the charges of the extremal BPS 
black hole and the condition of zero time-like NUT charge. 



To see why this is true, expand the coset element ksps that generates the BPS attractor 
flow using a aA . 

k B ps = a aA C aA (5.2) 

where C aA are conserved along the flow. On the other hand, the conserved currents in the 
homogeneous space are constructed by projecting the one-form valued Lie algebra g~ l ■ dg 
onto k, which gives the vielbein in the symmetric space: 

g- l dgW = a aA V aA (5.3) 

where V aA is conserved: 

i ( V « A 4> a ) = (5.4) 



dr 

Therefore, the expansion coefficients of ksps are 

In terms of the vielbein, the supersymmetry condition that gives the BPS geodesies are 
written as [33] : 

V aA z a = (5.6) 

That is: 

V« A <p a z a = C aA z a = (5.7) 

Define z a = e^zp, 

C aA = C A z a (gg) 

Therefore, the coset element ksps is expanded by the coset algebra basis a aA as ksps = 
Oj aA C z a . 



Note that ksps has five parameters (C A , z) where A = 1, . . . , 4. As will be shown later, 
z can actually be determined by (C A ) and moduli at infinity. So the geodesic generated 
by ksps is indeed a four-parameter family. It is easy to show that ksps is null, but more 
importantly, it is nilpotent: 

4p5 = (5.9) 
As will be shown later, ksps indeed gives the correct BPS attractor flow. 



27 



5.1.2 Constructing k NonB ps 



To construct the non-BPS attractor flow, one needs to find an element in the coset algebra 
distinct from kpps that satisfies: 

^NonBPS = (5.10) 

The hint again comes from the BPS generator. Note that kpps — a-aAP A z a can be written 
as: 

k B ps = e~ zL ^k BPS e zL ^ (5.11) 
where k BPS spans only the right four coset generators a 1A : 

k° BPS = a 1A C A (5.12) 

That is, kpps is generated by starting with the element spanning the four generators anni- 
hilated by the horizontal SL(2) raising operator L~^, then conjugating with the horizontal 
SL(2) lowering operator L^. And it is very easy to show that (k BPS ) 3 = which proves 
(k B ps) 3 0. 

In G 2 (2)/SL(2, R) 2 , there are two third-degree nilpotent generators in total. And since 
there are only two SX(2,R)'s inside H, a natural guess for a non-BPS solution is to look at 
vectors with fixed properties under the second SL(2, R) group. An interesting condition is to 
have positive charge under some rotation of L 3 , i.e. an SL(2, R) rotation of J2a=i 2 a *AC aA . 
Therefore, this suggests to us to start with the element spanning the four generators annihi- 
lated by the square of the vertical SL(2, R) raising operator (L+) 2 and then conjugate with 
the vertical SX(2,R) lowering operator L~: 

k N onBPs(z) = e- zL "k° NonBPS e zL " (5.13) 

where 

k NonBPS = a *aC aa where a, a = 1,2. (5.14) 

And one can show that: (k% onBPS ) 3 = which proves (kNonBPs) 3 = 0. Moreover, (kNonBPs) 2 
is rank two. 

As long as one can pick the coefficients C aA and the twistor z that describes the 5*0(1, 2) 
direction to be such that the time NUT charge is zero, this generator will give nice non-BPS 
extremal black holes. All the known non-BPS solutions may be recovered this way, and 
more, as this construction gives absolute freedom to pick the charges and moduli at infinity 
for the black hole (clearly for certain values of charges and moduli the solution will crash 
into a naked singularity, but this is to be expected from comparison with the BPS case) 
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5.2 Properties of Flow Generators 



5.2.1 Properties of kpps 

We now turn to solving for v a and w a in (14.401) in terms of C A and z. First, from (14.391) 
we know that the null space of k 2 is five-dimensional and the v a span the two-dimensional 
complement of this null space. For kpps = CLaAC A z a the null space of (kpps) 2 does not 
depend on C A . Therefore, the v a depend only on the twistor z = z 2 /z 1 . 



Recall that we are using the basis where k has the form ( 1 4.281) . From inspection of k BPS 
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we find that {vx,v 2 ) can always be chosen to have the formi 

vi = (V 1} -ViVt, 0) v 2 = (-V 2 , Vl V 2 , y/2) (5.17) 
where T]i is a 3d metric of signature (1,-1,-1), and V\, V 2 are two three- vectors with 

V 1 -V 1 = V 1 -V 2 = Q V 2 -V 2 = -l (5.18) 

Since any linear combination of (v\,v 2 ) forms a new set of (v\,v 2 ), this means in particular 
that any v 2 + cv\ gives a new v 2 . Looking at the forms of (vi,v 2 ), we see that V 2 is defined 
up to a shifting of V\ as V 2 = V 2 — cV\. 

An explicit computation shows that V\ and V 2 are given by the twistor z and u as 

/ (z l ) 2 + {z 2 f\ fz 1 u 1 + z 2 u 2 \ 

V x =\ (z') 2 -(z 2 ) 2 V 2 = l 2 L 2 J 2 V-A 2 . (5.19) 

V 2zh 2 ) z u - y z , u2 + z2ul J 

2 

where the twistor u = \ is related to c by 

1 + 2cz . 
u = -—z (5.20) 

The twistor representation of V\ and V 2 are 

Vf = 2z a z li Vf = z a vP + A" (5.24) 



4 When solving for (v,w), there are some freedom on the choice of (vi, V2) and (wi, W2): firstly, a rotational 
freedom 

(vi,v 2 ) (vi,v 2 ) ( p 12 ) and (w x , w 2 ) (w u w 2 ) ( J 1 p 12 ) (5.15) 

y K 2 i K22 J \ n 2 \ K 2 2 ) 

where R is orthogonal: RR T = 1 Secondly, a rescaling freedom: 

v a — > rv a and w a — * -w a (5.16) 

r 

5 With the inner product of three-vectors defined as 

Va-Vb = vlvi v b (5-21) 
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where we have used the rescaling freedom to set z x u 2 — z 2 u l to be 1. Note that for the BPS 
case, the twistor u is totally arbitrary. 

Now we solve for w a . The condition that w a are orthogonal to v a dictates that they have 
the form: 

wi = (W 1: mWu 0) w 2 = (W 2 , ViW 2 , 0) (5.25) 
where W\ and W 2 are linearly independent, and are related to the charges by w a ■ Wb = c a b- 

W x -W x = leu W 1 -W 2 = X -c X2 W 2 -W 2 = X -c 22 (5.26) 

Recall that V 2 is defined up to a shift by V\. V 2 = V 2 ° — cV\. The consequence is that W 1 is 
defined up to a shift by W 2 : W\ = W® + cW 2 . Note that the numerical factors in front of V\ 
and W 2 are opposite. Write down (W 7 ^, W 2 ) i n terms of (C A ,z): 

( (C 2 + C 4 ) + (C 1 + C 3 )z \ ( -(C 2 + C 4 ) + (C 1 + C 3 )z 

W? = j- [C 2 - C 4 ) + (C l - C 3 )z W 2 = -\ -{C 2 - C 4 ) + (C 1 - C 3 )z 

Az \ 2C 3 + 2C 2 z J 2 \ -2C 3 + 2C 2 z 

The twistor representations of W\ and W 2 are 



(5.27) 



CW-CPu 1 CV-CV\ _ / CV-CV CV-CV . 

"1 ~~ ( /^2„,2 r^3„,l r~<3„,2 /nr4„,l ^2 — I ^2 ^2 /ntfU ^3^2 ,^4 1 (5.ZOj 



CV-CV CV-CV j CV-CV CV-CV 

Define the totally symmetric P a ^\ 

P ni = c i p^ = C 2 P 122 = C 3 P 222 = C 4 (5.29) 
Then the three- vectors (W±, W 2 ) span the four dimensional space 

(W?, W«) BPS = (P a ^, P^z,) (5.30) 



5.2.2 Properties of k NonB ps 



The form of v a for the non-BPS case is only slightly different from the BPS case: the two 
vectors v a can be chosen to have the form: 

vi = (V u rjiVi, 0) v 2 = (V 2 , - m V 2 , V2) (5.31) 
The twistor representation of a three-vector v = (x, y, z) is 

a v = xa + ya 3 + zai = ^ X * V ^ Z _ y ^ (5.22) 

It's length is 

v T mv = det(a v ) = x 2 -y 2 - z 2 (5.23) 
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where Vi, V2 are two three- vectors satisfying the same condition as the BPS ones f !5. 18|) . 
Again, the vectors V\ and V 2 can be written as (I5.19p . and the twistor representations are 
given in ( 15.241) with one major difference: u is no longer arbitrary, but is determined by C aA 
as: 

u 2 C 22 , . 

u = ^ = ( 5 - 32 ) 

The form of {w\,w 2 ) are also slightly different from the BPS one (15.251) 

wi = (Wi, -ViW h 0) w 2 = (W 2 , V1W2, 0) (5.33) 
The (Wi, W 2 ) can be written in terms of (C aa , z) thus: 



U2\21 



1 



[(C U C 22 - C 12 C 2V )z 2 + (C 22 ) 2 ] + [C U C 22 - C 12 C 21 + (C 1 
W i = 0^22 r 7n2 2^ ( -l(C u C 22 - C l2 C 2l )z 2 + (C 22 ) 2 } + \C ll C 22 - C 12 C 21 + (C 12 ) 2 } 
2(0 z — C z)y 2[(C n C 22 -C 12 C 21 )z + C 12 C 22 } 

(5.34) 

z[C n z 2 + (3C 12 - C 21 )z - 2C 22 ] + [C u z + C 12 - C 21 ] \ 
W 2 = --\ -z[C u z 2 + (3C 12 - C 2l )z - 2C 22 ] + [C u z + C 12 - C 21 ] (5.35) 

2[C n z 2 + (2C 12 - C 21 )z - C 22 } ) 

In terms of u = ^ = §17, the twistor representation of W\ and W2 are: 

Wf = u a u? + (C u u 2 - C l2 u l )z a zP (5.36) 

= ( z a u p + u a z p ) + (C 21 - C 11 ^ - Sm 1 )^ (5.37) 

As a consequence, the precise value of u is an extra constraint on the vectors w a , and there 
is only a three-dimensional space of them, with (Wi,^) a linear combination of (0, Vi), 
(^,0) and (u a u^2V 2 ). 

(Wi, W 2 ) NonBP s = m(0, 7i) + 0) + t{u a vP, 2V 2 ) (5.38) 



6 Single-centered Attractor Flows in (2 2 (2)/(»SX(2, R) x 

SL(2,R)) Model 

Now that we have completely characterized the generators of single-centered attractor flow, 
we can lift the geodesies to four-dimensional black hole solutions. After some calculational 
effort, we find that the BPS solution is given in terms of harmonic functions. Next, we show 
that the non-BPS case is qualitatively different, and the final solutions cannot be formulated 
so simply. 
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6.1 BPS Attractor Flow 



6.1.1 Lifting a geodesic to AD 



We have already noted that the flow starting from generic asymptotic moduli (x ,y ) is 
generated by M(r) = e^ kr+9 ^ 2 , with g defined before in (14 .43 p . The matrix g has the same 
form as k. Therefore, in the BPS case, it has the expansion 



9bps — a a AZ a G 



(6.1) 



where the twistor z is the same as the one in k B ps- The flow of (x, y) and u can be extracted 
from the symmetric matrix S(t) = M(t)SqM(t) t via (14.321) and (I4.33p . Using k BPS = 
and g 3 = 0, S(t) = S e kr+9 is a quadratic function of r: 



S55( t )bPS 

S35 (t) bps 
S33 (t)bps 



a B (r) + (3 b {t) 
1b{t) + S B (r) 
£b(t)+( b (t) 



(6.2) 



where {a B (r), j b (t), e B (r)} are quadratic functions of r, and {{3b,8b,(b} are linear func- 
tions of r: 



a B {r 
Pb(t 

1b(t 
5b(t 

Cb(t 



z\H 2 H* - (H ) ) + z(H H - H H ) + (H H - (H 2 ) 2 ) 
(H 2 - H A )z + (H 3 - H 1 ) 

-- ((H 1 ^ - H 2 H 3 ){z 2 - 1) + 2{H 2 {H 2 + H A ) - H^H 1 + H 3 ))z) 
2 



--((H 1 + H 3 )z - (H 2 + tf 4 )) 



(H'H 3 - (H 2 ) 2 )z 2 + (H'H 4 - H 2 H 3 )z + (H 2 H 4 - (H 3 ) 2 ) 



2tt3\ 



^3^2^ 



= 2{H 2 z + H 3 ) 

where H A is a linear function of r defined as H A (r ) = C a t + G A . 



(6.3) 



The attractor values are reached when r — > 00 along the geodesic: 



x 



(k 2 S ) 



BPS 



(k 2 S ) 33 



Vbps 



V(k 2 s ) 33 (k 2 s ) 55 -(k 2 s ; 



2 

35 



(k 2 S { 



0J33 



and 



1 



u 



BPS 



) 2 



(6.4) 



(6.5) 



V(k 2 S ) 33 (k 2 S ) 55 - (k 2 S ) 

The asymptotic moduli (20,2/0) can be expressed in terms of (G A , z) by extraction from 
S = e 9 S . 
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Figure 3: Sample BPS flow. The attractor point is labeled (x*,y*). The initial points of 
each flow are given by (xi = 1.5, y\ = 0.5), (x 2 = 2, y 2 = 4), (x 3 = — 0.2, y 3 = 0.1), (x^ = 
—1,2/4 = 3) 
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Using this technique, one can construct all BPS single-centered black holes. The charges 
of each black hole can be read off from the current J using (I4.34p . One example is given 
in Figure [3], where we parametrically plot x(t) and y(r) for a BPS black hole with charges 
(p°,p\qi,q ) = (5,2,7,-3) and attractor point (x*,y*) = (0.329787,0.788503). E 



6.1.2 AD solution for given set of charges 

To get the solution for a specific set of charges requires more effort. In this section, we 
present the analytical result for any set of charges (p 1 , qi). 

The ten parameters in kpps and g are {z, u, C A , G A }, among which the twistor u is 
arbitrary, corresponding to the freedom of the shift by (fi,^) in the definition of (v 2 ,Wi): 
(t> 2 , Wi) — > (t>2 + cv i, W\ — cw 2 ). The remaining true parameters are enough to parameterize 
the four D-brane charges (p 1 , qi) and the arbitrary asymptotic moduli (x , y ) under the 
condition of vanishing Taub-NUT charge and fixing u = 1. We now solve for ksps and 
g for the given D-brane charges and (xo,yo), using the eight constraints, namely, 4 charges 
and zero Taub-NUT charge plus 3 asymptotic moduli, to fix C A and G A , leaving the other 
twistor z unfixed 

For the sake of simplicity, we will denote kpps by k for the rest of this section. Then the 

7 T 

current J(Q) = — gives the five coupled equations: 

Q = S (k + ±[k,g])S (6.6) 

In order to show that the BPS flow can be expressed in terms of harmonic functions: 

H(r) = Qr + h with Q = ( p I ,q I ) and h=(h I ,h I ) (6.7) 

we will solve g in terms of h instead of (xq, y^). The four /z's relate to the asymptotic moduli 
by 

x = x*(Q^h) y = y *(Q^h) u = u*(Q^h) (6.8) 
and there is one extra degree of freedom to be fixed later. 

To evaluate [k,g], we first use the commutation relation (14.111) to obtain 

[a 1A C A , a 1B G B ] = (C, G) (-4L+) (6.9) 

6 The discriminant of the charge (5, 2, 7, —3) is positive, so this is indeed a BPS solution. 
7 The twistor z can be left unfixed because we will not specify the asymptotic values of the scalars with 
translational invariance, namely, ({a, A 1 , Bi}). Fixing them can fix the twistor z. 
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where the product between C A and G A is defined as (C, G) = C l G 4 -3C 2 G 3 +3C 3 G 2 -C A G\ 
Then twisting Eq (16.91) with the twistor z as in (15.111) gives the commutator of k and g with 
the same twistor z: 

[k, g] = [a aA z a C A , a pB z p G B ] = (C, G)S (6.10) 
where 6 is defined as 6 = —^-^e~ zL ^L\e zL h . On the other hand, using (14.431) . 

[k, g] = (v 2 vj - v 1 v^)S (w 2 -rm-Wf m 2 ) (6.11) 

G can also be written as = (v 2 vf — ViV 2 )S , and we can check that {w 2 ■ mi — W\ ■ m 2 ) = 
(C,G). 

First, separate from G A the piece which has the same dependence on [h 1 z) as C A on 
(Q,z): 

G A = G A + E A with G A = C A (Q -> h,z) (6.12) 
That is, g contains two pieces: 

9 — 9h + A with g h = a a AZ a G A and A = a a AZ a E (6.13) 

We need to solve for E A . 

There are three constraints from (16.81) . The (xo>2/o) an< ^ u o are ex t rac ted from the sym- 
metric matrix S = e 9 So via (I4.32p and (14.331) . On the other hand, requiring (16.81) gives 
(x ,y ,u ) in terms of h: 

(glS ) 35 _ a/(^5 , o)3 3 (^5'o)55 - (glS )l 5 1 

'"' (^o) 33 m ~ (glSohs ^(^0)33(^0)55-^0)^ 

(6.14) 

Therefore, defining II = (e 9 — 2 )So, n has to satisfy three constaints: 

n 33 = 1135 = n 55 = (6.15) 

in order for (16.141) to hold for arbitrary h. Using the unfixed degree of freedom in h's to set 
(C, G h ) = 0, dSHJ) becomes 

Q = S (k + ^[k,A])S (6.16) 

The zero Taub-NUT charge condition in (16.161) imposes the fourth constraint on A: the 
(3,2)-element of So(k + %[k, A]) So for arbitrary k has to vanish. Combining with (16.151) . we 
have 4 constraints to fix E A to be: 

E l = -E 3 = -~^— 2 E* = -E* = -j— 2 (6.17) 
1 + z z 1 + z A 



35 



The remaining 4 conditions in the coupled equations (16.161) determine C A in the BPS 
generator ksps = &aAZ a C A to be 



c 1 


= V2- 


c 2 


= V2- 


c 3 


= V2 


c 4 





-% — (l\ z — Sp^-z 2 + p°z 3 



{1 + z 2 ) 2 

lp l - v )z + (p° + 2^)z 2 + p 1 z 3 

>° + 2f )z + (2p x - Vo )z 2 - f z 3 
(1 + z 2 ) 



2\2 



(6.18) 



The are then determined by G A = a a AZ A C A (Q — > fa, 2;). Using the solution of C A and 
G A , we see the product (C A , G A ) is proportional to the symplectic product of (p 1 , 5/) and 
(/''•/'/): 



(C" 4 , G^ 4 ) = - — — ^ < Q,h > where < Q, fa >= p°h + p x h x - q x h x - q h° (6.19) 
1 + z 2 

The condition (C A , G A ) = is then the integrability condition on fa: 

<Q,h>= p°h + p l h x - qih 1 - q h° = (6.20) 



Substituting the expressions of C A and G A in terms of (p^qi) into (16. 2p . we obtain the 
BPS attractor flow in terms of the charges (jp 1 ,qi)- In particular, the attractor values are 



2[(^+p°f] 2[(pT+P°f] (b - 21) 

where J^p ,]? 1 , q\, qo) is the quartic £7(7) invariant: 

J4(p°,p\qi,qo) = 3(Ai) 2 - 6(p%)(Ai) - (p%) 2 - 4(p 1 ) 3 9o + 4p°( gi ) 3 (6.22) 

thus J^p^p 1 , y,<?o) is the discriminant of charge. The attractor values match those from 
the compactification of Type II string theory on diagonal T 6 , with q 1 — > y. The attractor 
value of w is 

«*bps = (6.23) 

The constraint on h from u = 1 is then J^(h°, fa 1 , — , fao) = 1- 

Now we will show that the geodesic we constructed above indeed reproduces the attractor 
flow given by replacing charges by the corresponding harmonic functions in the attractor 
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moduli. Using the properties of A, we have proved that, in terms of k and g, the flow of 
(x, y) can be generated from the attractor value by replacing k with the harmonic function 
kr + g h : 

x (r)= x *(k^kr + g h ) y{r) = y*{k -> kr + g h ) (6.24) 

Since k and gh have the same twistor z, this is equivalent to replacing the C A with the 
harmonic function C A r + G A while leaving the twistor z fixed: 

x (r) = X *(C A - C A r + G A , z) y{r) = y*(C A -> C A r + Gjf, 2) (6.25) 

Since G A is linear in Q and G A linear in h, and since z drops out after plugging in the solution 
of C A in terms of (Q, z) and G A in terms of (h, z), this proves that the flow of (x , y ) is given 
by replacing the charges in the attractor moduli by the corresponding harmonic functions: 

Xbps(t) = x* B ps(Q ^Qr + h) Vbps(t) = Vbps(Q ~* Q T + h ) ( 6 - 26 ) 
The integrability condition < Q,h >= 0, in terms of H = Qt + h, is 

<H,dH>=0 (6.27) 



6.2 non-BPS Attractor Flow 
6.2.1 Lifting a geodesic to AD 

Similar to the BPS attractor flow, the non-BPS flow is generated by M(r) = e^ kr+9 ^ 2 , 
and (x, y) can be extracted from the symmetric matrix S(t) by (I4.32p and the relevant 
elements of S(t) are given by ( 16. 2ft . The only difference is that now {as, Pb, 1b, 8b, £b, (b} 
are changed into the non-BPS counterparts {ctNB, Pnb, Inb, &nb, £nb, (nb}, which can be 
written in terms of H aa = C aa r + G aa and z: 

a NB (r) = -((H 21 H 12 -H 22 H u )(z 2 -l) 2 + (H 22 ) 2 z 2 -2zH 12 H 22 + (H 12 ) 2 ) 

(3 NB { T ) = (z 2 -l)(H 11 -H 21 z)-3z 2 H 22 + 2zH 12 + H 22 

lNB (r) = (z 2 -l)(2z(H u H 22 -H 12 H 21 ) + H 22 H 12 ) + z((H 22 ) 2 -(H 12 ) 2 ) 

5 nb (t) = ^(z(l + z 2 )H 11 + z 2 (3H l2 -H 21 )-2H 22 z+(H l2 -H 21 )) 

e NB (r) = (AH u H 22 -H 12 (H 12 + AH 21 ))z 2 -2H l2 H 22 z-(H 22 ) 2 

( NB ( T ) = 2(H u z 2 + (2H 12 + H 21 )z + H 22 ) (6.28) 

Note that = u is fixed, independent of r. The non-BPS flow written in terms of (H aa , z) 
has the same simple form as the BPS flow, i.e. the scalars are rational functions with both 
the numerator and denominator being only quadratic. This is due to the nilpotency of the 
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generator: k 3 = 0. Again, the attractor values are reached when r — > oo, and the asymptotic 
moduli can be expressed in terms of (G aa , z) by extraction from S = e 9 So. 

Unlike the BPS case, there are only eight parameters in UnotiBps and QnouBps'- the two 
twistors (z,u) and (C aa ,G aa ) under the constraints that 



/J22 Q22 



Therefore, while ksps and gsps can parameterize arbitrary (p 1 , qj) and (x , yo) while leaving 
(z,u) free, all the parameters in k^ on BPs and gNonBPS, including (z,u), will be fixed. 




Figure 4: Sample non-BPS flow. The attractor point is labeled (x*,y*). The initial points of 
each flow are given by: (xi = 0.539624, yi = 5.461135), (x 2 = 1.67984, y 2 = 0.518725), (x 3 = 
-0.432811,2/3 = 0.289493), (x 4 = 1.28447, y 4 = 1.49815), (x 5 = -0.499491, y 5 = 0.181744) 



The attractor point in terms of C aa is 

(6.30) 
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with z given by 

I (_3C 12 - C 21 + 2 {{z 2 - 3)C U + 3^(C 12 + C 21 ) + 6C 22 )) = (6.31) 

As in the BPS case, the charges of the black hole are read off from the current J using ( 14.34[) . 
We have checked that the attractor point is a non-supersymmetric critical point of the black 
hole potential V BH = \Z\ 2 + \DZ\ 2 : 

dV BH = and DZ ± (6.32) 

It reproduces the results reported in the literature [9]. An example of the non-BPS attractor 
flow is shown in Figure HI with (p°,p 1 ,qi,qo) = (5,2,7,3) and attractor point (x*,y*) = 
(-0.323385, 0.580375). Note that J 4 (5, 2, 7/3, 3) < 0, so this is indeed a non-BPS black hole. 
Unlike the BPS attractor flow, all the non-BPS flows starting from different asymptotic 
moduli have the same tangent direction at the attractor point. The mass matrix of the black- 
hole potential at a BPS critical point has two identical eigenvalues, whereas the eigenvalues 
at a non-BPS critical point are different. The common tangent direction for the non-BPS 
flows corresponds to the eigenvector associated with the smaller mass. 



6.2.2 AD solution for given set of charges 



We now discuss how to construct the non-BPS black hole solution for a specific set of charges 

(p^qi)- 

One major difference between the non-BPS case and the BPS case is that 

[kNonBPS, dNonBPs} = (6.33) 

automatically, since the forms of (u>i, io 2 ) and (mi, m 2 ) guarantee that w± ■ m 2 = w 2 ■ ra x = 0. 
Thus the charge equation (16.61) becomes simply 

QnouBps = So(kN on BPs)So (6.34) 

These five coupled equations determine the two twistors (z,u) and C aa in terms of (p^qi)- 
Similar to the BPS case, the four equations which determine the D-brane charge allow us to 
write C aa in terms of the charges (p ,^ 1 , qi, qo) and the twistor z via 

(-2q + 6(p x - q )z 2 + 4(p° + gi)^ 3 - 6pV) 



C 21 



v / 2(l + ^ 2 ) 3 

(p° + f ) - 2(2p x - q )z - (p° + 5f )z 2 + 2p 1 z 3 
v^l + z 2 ) 2 

(p° - qi) - ^q\Z 2 + 4(3p 1 - q )z 3 + (3p° + q^z 4 



\/2(l + 



,2\3 



c22 2p l + (pP + 5f)z - 2(2p 1 - q )z 2 - (p° + | )z* 
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and u = ^jtj. In contrast to the BPS case, the G aa do not enter the equations and therefore 
cannot be used to eliminate the twistor z. Requiring the Taub-NUT charge to vanish gives 
the following degree-six equation for the z: 

p°z 6 + Qp l z^ - (3p° + 4gi)/ - 4(3p 1 - 2q )z 3 + (3p° + A qi )z 2 + 6p 1 z - p° = (6.36) 

The three parameters in gNmBPS-i namely, G aa with the constraint ^ = u are then fixed 
by the given asymptotic moduli (x , yo) and u = 1. 

Similar to the BPS flow, the full non-BPS flow can be generated from the attractor 
value by replacing C aa with the harmonic function H aa (r) = C aa T + G aa , while leaving 
z unchanged as in (16.251) . However, there are two major differences. First, the harmonic 
functions H aa have to satisfy the constraint 

H 2 \r) _ _C^_G^_ 

Note that this does not impose any constraint on the allowed asymptotic moduli since there 
are still three degrees of freedom in G aa to account for (xq, yo, «o)- We will see later that it 
instead imposes a stringent constraint on the allowed D-brane charges in the multi-centered 
non-BPS solution. 



Secondly, unlike the BPS flow, replacing C aa in the attractor moduli by the harmonic 
function H aa (r) is not equivalent to replacing the charges Q with H = Qr + h as in (I6.26p . 
The twistor z here is no longer free, but is determined in terms of the charges as a root of the 
degree-six equation (16.361) . so replacing Q by Qr + h, for generic Q and h, would not leave 
z invariant. Therefore, the generic non-BPS flow cannot be given by the naive harmonic 
function procedure, as proposed by Kallosh et al [19J. Next, we will define the subset of the 
NonBPS single-centered flow that can be constructed by the harmonic function procedure. 

When the attractor has only DA — DO charges, namely, = (0,p 1 ,0,qo), (I6.36j) has a 
root z = 0, which is independent of the value of charges. If the asymptotic moduli h is also 
of the form of ^ 40 = (0, h 1 , 0, h ), replacing Q 40 by Q^t + /i 40 would leave the solution z = 
invariant. Now we will use the duality symmetry to extend the subset to a generic charge 
system with restricted asymptotic moduli. 

The one-modulus system can be considered as the STU attractor with the three moduli 
(S,T,U) identified. Since the STU model has SX(2,Z) 3 duality symmetry at the level of 
the equations of motion, the one- modulus system has an SL(2, Z) duality symmetry coming 
from identifying the three SL(2, Z) symmetries of the STU model. That is, the one-modulus 
system is invariant under the following element of SL(2, Z) 3 
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The modulus z = x + iy transforms as 



Tz 



az + b 
cz + d 



(6.39) 



and the transformation on the charges is given by [35]. A generic charge (p°,p 1 ,qi,qo) can 
be reached by applying the transformation T on a DA — DO system. 

Under the aforementioned transformation, a DA — DO system transforms into f Q40 

/ \ 



Q 



10 



p 




TQ 



40 



V ft / 

The solution of the twistor z with the new charges VQm is 



/ -c(3dY + c 2 q ) \ 
d(2bc + ad)p 1 + ac 2 qo 
3(b(bc + 2ad)p 1 + a 2 cq Q ) 
a(Sb 2 p l + a 2 q ) 



± Va 2 + c 2 



(6.40) 



independent of the D4 — DO charges we started with. Now given an arbitrary charge Q, 
there exists a transformation Tq such that Q = TqQ^ for some Q^q. The twistor z remains 
invariant under Q — ► Qr + Tq^o for arbitrary /i 40 . We conclude that the non-BPS single- 
centered black holes that can be constructed via the naive harmonic function procedure are 
only those with (Q, h) being the image of a single transformation T on the (Q40, ^40) from a 
DA - DO system: 



X NB {T 



X 



NB 



(f q 40 — > f <5 4 o^ + r/i 



10 J 



Z/jvb(^) = Z/vs(rQ4o -> f Q 40 r + Th 40 ) (6.41) 



Since we are considering arbitrary charge system, the constraint is on the allowed values of 
h. 



7 Multi-centered Attractor Flows in G 2 ( 2 )/{SL(2,R) x 

SX(2,K)) model 

As proven in the pure gravity system, the multi-centered attractor solutions are given by 
exponentiating the matrix harmonic function K(x)\ 

S(x) = e K{s) S (7.1) 

with K(x) having the same properties as the generator k: 

K 3 (x) = and K 2 (x) rank two (7.2) 
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We now describe how to formulate K(x) for multi-centered solutions in G 2 (2)- 

The K(x) satisfying all the above constraints is constructed as 

K(x) = KW a (x) T + W a (x) V T]S (7.3) 

0=1,2 

with v a being the same two constant null vectors in k, and the multi-centered harmonic 
function 

W a (x) = T-^- + m a (7.4) 

i ' * ' 

is everywhere orthogonal to v a . The two 7- vectors (m 1 ,m 2 ) contain the information of 
asymptotic moduli and has the same form as (wi, w 2 ). Write K(x) as K(x) = J2i ig-g-l 9 
where 

k i = ^2l v a( w a)J + (w a )ivJ]S and g = Y,[v a ml + m a v£\S (7.5) 

0=1,2 0=1,2 

Since v only depends on the twistor (z,u), and (wi, w 2 ) are linear in C A or C aa , the above 
generating procedure is equivalent to replace C A or C aa by the multi-centered harmonic 
functions while keeping the twistor (z, u) fixed. 

Next we discuss the properties of the BPS multi-centered attractor solution and non-BPS 
ones separately, since they are very different in character. 



7.1 BPS Multi-centered Solutions 



In constrast with the multi-centered solutions in pure gravity, now the second term of the 
current J = VK + ~ [VK, K] does not vanish automatically since 

Therefore, the centers are no longer free, and we cannot simply read off the charges from J. 
Instead, we need to solve for C A and G A in a set of 5iV coupled equations. The divergence 
of the current is 

V • J = 4tt 8(2 ~ Zi)So(ki + ~ [K 9] + lY F^4) 5 ° ^ 

i 

Using Qj to denote the charge matrix which relates to the D-brane charge {p ^ 1 , qi, qo}i as 
in (14.361) . and with Q32 as the vanishing NUT charge, we have 5N coupled equations from 

Qi = S (ki + hi*, g] + \Y J^L)5 (7.8) 

3 J 
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The generators of the multi-centered BPS attractor solution {ki} and g have 4(jV + 1) + 2 
parameters in total: the two twistors (z,u) and {Cf,G }. On the left hand side of (17. 8p . 
there are also 3iV — 3 degrees of freedom from the position of the centers Xj. On the other 
hand, a generic iV-centered attractor solution has 4N D-brane charges (p^qi), and three 
additional constraints from the asymptotic moduli (xo,yo) and uq = 1. As we will show, like 
the single-centered BPS solution, the three asymptotic moduli, together with the vanishing 
of the total Taub-NUT charge, determine the 4 G A inside g. Moreover, as in the single- 
centered case, we can solve Cf in terms of the 4D D-brane charges Qi while leaving (z, u) 
unfixed. The remaining N — 1 zero Taub-NUT charge conditions at each center will impose 
N — 1 constraints on the distances between the iV-centers Xi. 

First, integrating over the circle at the infinity, £\ Qj = j- J V ■ J gives the sum of the 
above N matrix equations: 

Qtot = J2^ = 5 »(E ki + 5 E 3\) S o (7-9) 

i i i 

which is the same as the one for the single-center attractor with charge Qtot- This determines 
g to be g = gh + A, same as the one for single-centered attractor in (16.1 3D . As in the single 
centered case, h is fixed by the asymptotic moduli (xq, yo) by 

x o = x*bps(Q vo = v*bps(Q h ) ( 7 - 10 ) 

and the two constraints: 

<Q toU h>=0 J 4 (h°,h\ f j,h ) = l (7.11) 

We have used the vanishing of the total Taub-NUT charge to determine A. Next, we will 
use the remaining coupled 5A^ — 1 equations to solve for the AN {Cf~} and impose — 1 
constraints on the relative positions between the N centers {x^} where i = 1, • • • , N. 

The tentative solutions of Q 4 are given by (16.181) with {p 1 , qj) replaced by (pj , qi,i). The 
flow generator of each center ki is then hi = a a AZ a C{ L . Substituting the solution of k{ and 
g = gh + A into (17.81) , and using 

[ki, g h ] = 2<Q l ,h>Q [h, kj\ =2<Qi, Q j > 6 (7.12) 

where all the fcj's and gt have the same value for the twistor z, we get 

Qi = 3 (ki+ <Q l ,h>® + hki, A] + J2 ^%^ Q)^o (7.13) 

j 

Just as in the single-centered case, the solution of ki and the form of A guarantee that 

Q l = S (k l + ^[k i ,A})S (7.14) 
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We see that as long as the following integrability condition is satisfied: 



<Q i ,h>+Y) < l T Q >> = (7.15) 

^^^^ / >-» . rp . \ 

3 J 

the ki and g given above indeed produce the correct multi-centered attractor solution. Just 
like in the single-centered case, the multi-centered solution flows to the correct attractor 
moduli (x*,y*) near each center, independent of the value of z. It also follows that the 
multi-centered solution can be generated by replacing the charges inside the attractor value 
by the multi-centered harmonic function: 

xbps(x) = x* BPS (Q -> Y] i + h) Vbps(x) = y* BPS {Q -> Y] rJ^r, + h ) ( 7 - 16 ) 



The sum of the iV equations in the integrability condition (17.151) reproduces the constraint 
on h: < Qtot, h >= 0. Thus the remaining N — 1 equations impose N — 1 constraints on the 
relative positions between the N centers {x*i} with % — 1, • • • , N. From (1A.5I) and (14.61) . we 
see that *duj is given by J23. Defining the angular momentum J by 

uji = 2e ijk J j — as r — >• 00 (7.17) 
we see that there exists a nonzero angular momentum given by 

Kj 

Thus we have shown that our multi-centered BPS attractor solution reproduces the one 
found in US]. 



7.2 non-BPS Multi-centered Solutions. 



For given (z, u) and {C" a ,G act }, the non-BPS multi-centered solution is the same as the 
single-centered one as in (16.281) with H aa (T) replaced by the multi-centered harmonic function 
H aa (x) = J2i + G aa satisfying the constraint 

H} 2 {x) CP G 12 1 } 

Accordingly, the attractor values at each center is the same as (16.301) with the corresponding 
C" a . The asymptotic moduli are obtained by extraction from 5* = c 9 Sq. 
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The equation of motion for the non-BPS multi-centered solution simplifies a great deal 
since 

[VK(x),K(x)]=Q (7.20) 
automatically, following from the fact that for the non-BPS system: 

(wi)i ■ m 2 = (w 2 )i ■ mi = (7.21) 

which are guaranteed by the forms of NonBPS (wi, w^ji and (m l5 m 2 ). Therefore, the 5N 
equations (17. 8p decouple into iV sets of 5 coupled equations: 

Q, = S (k t )S (7.22) 



Equation (I7.22p differs greatly from the BPS counterpart (|7.8|) . Firstly, the generators of 
the multi-centered non-BPS attractor solution {fcj} and g have 3(iV+ 1) + 2 parameters: the 
two twistors (z,u) and {Cf a , G aa } with the constraint ( 17.191) . In constrast to the BPS case, 
g does not enter the equation. Thus we can simply use the three asymptotic moduli, without 
invoking the zero Taub-NUT condition, to determine the 3 G aa inside g. Secondly, unlike 
the BPS multi-centered solution, the position of the centers do not appear in the equation, 
therefore there will be no constraint imposed on them: the centers are free. Last but not least, 
the remaining 3iV + 2 parameters in (z, u) and C aa are not enough to parameterize a generic 
iV-centered attractor solution, which has 4N D-brane charges (pl,qi,i)- Accordingly, the 
multi-centered non-BPS attractor generated by this ansatz will not have arbitrary charges. 
Combining with the fact that Xi do not appear in the R.H.S of the equation, we find that 
all the iV vanishing Taub-NUT charge conditions can only act on the charges on the L.H.S. 
We conclude that, in total, there will be 2N — 2 constraints on the allowed charges. 

Now we will show in detail the derivation of the constraints. First, like in the single- 
centered NonBPS solution, the absence of the Taub-NUT charge at infinity fixes z via 

J2Qi = So(J2 k i) S ° ( 7 - 23 ) 

i i 

The solution is the same as the solution to ( 16.361) with the charges replaced by the total 
charges of N centers: z = z(Q — > YliQi)- Since all the N centers share the same twistor 
z, the absence of the NUT charge at each center imposes N — 1 constraints on the allowed 
charges Qf. all z(Qi) have to be equal. 

The remaining AN equations in (17.221) determine C aa in terms of z and Qi. Since the 
N-centers decouple, ( I7.22p for each center is the same as the single-center one (16.341) . Thus 
the solution of Cf a is given by (16.351) with (p 1 , qj) replaced by (pf, qi t i). Again, since all the 
centers share the same twistor u, the condition (17.191) imposes another N — 1 constraints on 
the allowed charges. Solving these 2N — 2 constraints, we see all the charges {Qi} are the 
image of a single transformation T on a multi-centered D4 — DO system Q^q^. 

Q l = f Q 40>4 (7.24) 
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The charges at different centers are all mutually local 

(Qi,Q,) = (7-25) 

Except for the constraint on the charges, the N centers are independent, and there is no 
constraint on the position of the centers. A related fact is that the the angular momentum 
is zero. 

Like the non-BPS single-centered case, though the multi-centered solution can be gener- 
ated from the attractor value by replacing C aa with the multi-centered harmonic function 
H aa {x) = +G aa under the constraint (I7.19p . while leaving z unchanged as in (I6.25p . 

the generic solution cannot be generated via the harmonic function procedure used in the 
BPS case, namely, by replacing the charges inside the attractor value by the corresponding 
multi-centered harmonic functions. The reason is again due to the fact that the twistor z, 
being a function of charges, does not remain invariant under this substitution of charges by 
harmonic functions. The multi-centered non-BPS solutions that can be generated by the 
harmonic function procedure are those with {Qi, h} being the image of a single T on the 
{Qw,i, ^40} of a pure D4 — DO system: 

xnb(x) = x* NB (TQ m -> |^ 40 :' + f /140) Vnb(x) = y* NB (TQ m -> V" . + f h m ) 

%Aj %aj ^ *A/ »X> £ 

i i 

(7.26) 

It appears that the existence of a simple linear ansatz for "superimposing" single center 
solutions exists in general only for mutually local extremal black holes, and only in the 
super symmetric case does it extend to mutually non-local centers. 

To summarize, the non-BPS multi-centered solution is different from the BPS case be- 
cause it imposes no constraints on the position of the centers, but instead on the allowed 
charges Qf. the choice of charges at each center are restricted to a three-dimensional sub- 
space, and they are mutually local. The result is that the centers can move freely, and there 
is no angular momentum in the system. It does not have interesting moduli spaces of centers 
with mutually non-local charges, so it is as "boring" as the pure gravity case. 



8 Conclusion and Discussion 



In this paper, we find exact single-centered and multi-centered black hole solutions in theories 
of gravity which have a symmetric 3D moduli space. The BPS and extremal non-BPS single- 
centered solutions correspond to certain geodesies in the moduli space. We construct these 
geodesies by exponentiating different types of nilpotent elements in the coset algebra. Using 
the Jordan form of these nilpotent elements, we are able to write them down in closed explicit 
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form. Furthermore, we can use a symmetric matrix parametrization to recover the metric 
and full flow of the scalars in four dimensions. 

We have also generalized the geodesies to find solutions for non-BPS and BPS multi- 
centered black holes. The BPS multi-centered solution reproduces the known solution of 
Bates and Denef. Given our assumption that the 3D spatial slice is flat, we find that a non- 
BPS multi-centered black hole is very different from its BPS counterpart. It is constrained 
to have mutually local charges at all of its centers and therefore carries no intrinsic angular 
momentum. It is possible that if we dropped this assumption, we could find more general 
non-BPS multi-centered solutions. Such configurations would probably be amenable to exact 
analysis only in the axially symmetric case, using inverse scattering methods. 

There are many other avenues for future work. One could explore nilpotent elements in 
other symmetric spaces, and see whether non-BPS bound states with nonlocal charges exist. 
In particular, it would be interesting to study E 8 e 8 \/ SO* (16), which is the 3d moduli space 
for d = 4, A/" = 8 supergravity. We would also like to find a way to modify our method 
so that we can apply it to non-symmetric homogeneous spaces, and eventually to generic 
moduli spaces. We could then study the much larger class of non-BPS extremal black holes 
in generic M = 2 supergravities. 
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A Derivation of the moduli space Ai% D 



Here we briefly review the derivation of the 3c? moduli space Ai^D from the c*-map of the 
AD supergravity coupled to ny vector multiplets [31] 132] [33] . 



The bosonic part of the action for the M = 2 supergravity coupled to ny vector-multiplets 



is: 



S 



16tt 



d i x^fg^ R - 2g i] dz i A * A dz j - F 1 A d 



(A.l) 



where the ranges of the indices are i,j = l,...,ny and / = 0,1,..., ny, and Gj 
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(ReM)ijF J + (ImJ\f)ij * F J . The complex symmetric matrix A/}j is defined by 

Fj = M LJ X J D t Fj = AT/jAX j (A.2) 

For model endowed with a prepotential F(X), 



(ImF ■ X) I (lmF ■ X)j /A . 

A/ " = f " + a x.£f.x (a - 3) 



where F u = didjF(X). 



After reduction on the time-like isometry, the action is S = — f dt f c? 3 x C The 3D 
lagrangian C has three parts: C = C gravity + C moduU + £ e . m where 



-gramtj/ = -^v/g R + d£7 A *d£7 - ^e 4U du A *dcu 
^mod«ii = ftjdz* A *dz J (A.4) 
£ e . m . = ^e- 2U (ImAf) IJ dA I A *dA^ + ^e 2l/ (/mAT)/j(dA / + A^duj) A *(dA J + A J du) 
+ (ReM) IJ dA I A (dA J + A J Q du) 

The dual scalars for 00 and A 7 are defined by: 

e 2U (ImM)jj * (dA J + A J dLu) + (ReM) u dA J = -d0 A , 

e 4C/ * dcu + (4d0 A / - <j> A idAl) = -d^ (A.5) 

After renaming the variables W — > a, Aq — > A 7 , A / — > £>/, we obtain the 3d lagrangian in 
terms of scalars only: 

£ = --^gR + dU A*dU + -e- 4U {da + A^Bj - BjdA 1 ) A*(da + A^Bj - BjdA 1 ) 
z)d^ A *dz j + -e~ 2U (ImJ\f)ijdA I A *dA J 
+ l -e- 2U {ImN~ l ) IJ {dB I + (iteA/VdA*) A *(dBj + {ReM) JL dA L ) 
= -^R + g mn d a <p m d a <p n (A.6) 

where, as before, <p n are the 4(ny + 1) moduli fields: n = {£/, z 1 , z\ a, A 1 , _B/}, and g ab is 
the space time metric, g mn is the moduli space metric. Therefore, the moduli space M.?,d 
has metric: 



ds 2 = dU dU+ ^e- 4U (da + A 7 d5j - B^A 1 ) ■ (da + A 7 dfi 7 - B/dA 7 ) + g fj (z, z)dz l ■ dz j 



+^e~ 2U [(ImAf- 1 ) IJ (dB I + N IK dA K ) ■ (dBj + JT JL dA L )] (A.7) 
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It is a para-quaternionic-Kahler manifold. Since the holonomy is reduced from SO {Any + 
4)) to Sp(2,M) x Sp(2n v + 2,R), the vielbein has two indices (a, A) transforming under 
Sp(2,M.) and Sp(2ny + 2, R), respectively. The para-quaternionic vielbein is the analytical 
continuation of the quaternionic vielbein computed in [31] : 



V 



a A 



III 

„a 



V \ 

iE a 



\ —v iu J 



The 1-forms are defined as 

u 
e a 



V = 



e K/2 - u X\dB I +N u dA J ) 
e a t dz l 

e- u e a i g f] e K/2 D 3 X I {dB I +M IJ dA J ) 
-dU + l -e~ 2U (da + A f dBj - B^A 1 ) 



(A.8) 



where e° is the veilbein of the 4D moduli space, and the bar denotes complex conjugate. 
The line element is related to the vielbein by 

ds 2 = -u-u + g al e a ■ e b - g al E a ■ E l + v ■ v = e a(3 e AB V aA ® V pB (A.9) 

where e Q/ g and €ab are the anti-symmetric tensors invariant under 5p(2, R) = SL(2,R) and 
Sp(2n v + 2,R). 

The isometries of the descends from the symmetry of the 4D system. The gauge 

symmetries in 4D gives the shifting isometries of -M^: 



A 1 — ► A 1 + AA 1 
Bj — > Bj + AB T 
a — > a + Act + AB I A I - AA 7 ^ 



(A.10) 



The conserved currents and charges are given by (14.51) and the discussion thereafter. 
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